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Abstract. We study integrability by quadrature of a spatially flat Friedmann model 
containing both a minimally coupled scalar field tp with an exponential potential 
V{ip) ~ exp[— y/6crK(p], k — ^8ttGn, of arbitrary sign and a perfect fluid with 
barotropic equation of state p = (1 — h)p. From the mathematical view point the 
model is pseudo-Euclidean Toda-like system with 2 degrees of freedom. We apply the 
methods developed in our previous papers, based on the Minkowsky-like geometry for 
2 characteristic vectors depending on the parameters a and h. In general case the 
problem is reduced to integrability of a second order ordinary differential equation 
known as the generalized Emden-Fowler equation, which was investigated by discrete- 
group methods. We present 4 classes of general solutions for the parameters obeying 
the following relations: A. a is arbitrary, h = 0; B. cr = 1 — h/2, < h < 2; CI. 
a = 1 - h/4, < h < 2; C2. a = |1 - h\, < h < 2, h ^ 1,4/3. We discuss the 
properties of the exact solutions near the initial singularity and at the final stage of 
evolution. 
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1. Introduction 

A large class of theories (multidimensional 0, |!|, Kaluza-Klein models, supergravity 
and superstring theories, higher order gravity, see for instance deal with a weakly 

coupled scalar field ip with a potential of the form 

V(<p) = ^e"^, (1) 



where k = \Z8ttGn, cr is a dimensionless positive constant, characterizing the steepness 
of the potential, and the constant Vo may be positive and negative. A number of 
authors ( see a l so references therein) have studied a spationally homogeneous 



and isotropic Friedmann model containing both a scalar field and a perfect fluid subject 
to the linear equation of state 

p = (1 - h)p, (2) 

where the constant h satisfied < h < 2. The attention was mainly focussed on the 
qualitative behaviour of solutions, stability of the exceptional solutions to curvature 
and shear perturbations and their possible applications within the known cosmological 
scenaria such as inflation and scaling ("tracking") . In particular, it was found by a 
phase plane analysis |[ f| that for "flat" positive potentials (Vo > 0, < a 2 < l — h/2) 
there exists an unique late-time attractor in the form of the scalar dominated solution. 
It is stable within homogeneous and isotropic models with non-zero spatial curvature 
with respect to spatial curvature perturbations for cr 2 < 1/3 and provides the power-law 
inflation. For "intermediate" positive potentials (Vo > 0, 1 — h/2 < a 2 < 1) an unique 
late-time attractor is the scaling solution, where the scalar field "mimics" the perfect 
fluid, adopting its equation of state. The energy- density of the scalar field scales with 
that of the perfect fluid. For h > 4/3 this solution is stable within generic Bianchi 
models to curvature and shear perturbations and provides the power-low inflation. The 
scaling solution is unstable to curvature perturbations, when < h < 4/3, although it 
is stable to shear perturbations. Regions on (cr 2 ,h) parametrical plane corresponding 
to various qualitative evolution are presented on figure 1. 

Integrability of the model is not well studied yet. Only for the special case for h = 1 
(dust) and a = 1/2 the procedure of getting a solution to the model was given in |TT 



In the present paper we study the problem of integrability by quadrature of a 
spatially flat Friedmann model containing both a minimally coupled scalar field with 
the exponential potential ([!]) and a perfect fluid with the equation of state (0) for several 
classes of sets of (cr, /i)-parameters. We apply the methods developed in our previous 



papers [12|-[I1| devoted to integrability of multidimensional cosmological models. It 
is clear that the possibility of reducing the problem to quadrature depends on the 
parameters a and h. We show that in a general case the problem is reduced to 
integrability of a second order ordinary differential equation known as the generalized 
Emden- Fowler equation, which was investigated by discrete-group methods HT5|,fTfi 



We present here 4 classes of general solutions for the parameters obeying the following 
relations. 
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Figure 1. The domains la, lb, Ha and lib are bounded by the lines h = 0, h = 
2, a 2 = 0, a 2 — 1/3, a 1 = 1 — ft/2. The solutions corresponding to la and Ila are 
inflationary on late times with an attractor stable to curvature perturbations. The 
late-time attractor for la and lb is the scalar dominated solution and for Ila and 
lib it is the scaling solution. The dotted curves A,B,C1 and C2 are identified in the 
text and present the integrable by quadrature classes of the model. 



A. a is arbitrary, h = 0. 

B. a = 1 - h/2, <h<2. 
CI. (7=1 — h/A, 0<h<2. 

C2. a = |1 - h\, < h < 2, h ^ 1,4/3. 
The corresponding curves on (h,cr 2 ) parametrical plane are presented on figure 1. 

The paper is organized as follows. In section 2 we describe the model and obtain 
the Einstein-scalar field equations in the form of the Lagrange-Euler equations following 
from some Lagrangian. Dynamical system described by the Lagrangian of this form 
belongs to the class of pseudo-Euclidean Toda-like systems. To integrate them we apply 
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in section 3 the methods developed in our papers [|T^, [13|, [L4| on multidimensional 
cosmology. The method used in the general case C (a 2 ^ (1 — h/2) 2 , < h < 2) 
is based on reducing the Euler-Lagrange equations to the generalized Emden-Fowler 
(second-order ordinary differential) equation. To conclude the paper, we discuss in 
section 4 the physical properties of the obtained exact solutions. Asymptotic behaviour 
of these solutions at early and late times is analyzed. Example of the explicit solution 
in the cosmic time is presented. 



2. The model 



We start with Einstein Eqs. in a spatially flat Friedmann metric 

ds 2 = _ e Mt) dt 2 + e M*) [ d r 2 + r 2 (d6 2 + sin 2 #d0 2 )] , (3) 

where exp[x(t)] = a(t) is the scale factor and the function a(t) determines a time gauge 
(a(t) = corresponds to the cosmic time t c gauge). We assume that the universe 
contains both a self-interacting scalar field p with a potential V(ip) and a separately 
conserved perfect fluid with the barotropic equation of state (0). The governing set of 
equations, which follows from coupled Einstein and scalar field equations, reads 

~ 2 a 



X 



X 



K 

y 



-if 2 + (V(p) + p) e 



2a 



la 



3^37 J X 



K 

T 



-^-{V{ip)-p)e 



2<i 



V{<p)e 



2a 



(p = (a — 3x)(p 
p = -3x(p + p). 

Using the barotropic equation of state we get immediately ^from Eq.(^) 



9 



PO -3(2-h)x 



(4) 

(5) 

(6) 
(7) 

(8) 



where po is an arbitrary positive constant. 

It can be easily verified, that the set of Eqs. (f|) -(H), where the presence of a perfect 
fluid density p and its pressure p is cancelled by Eqs.(|2])and (^), is equivalent to the set 
of Euler-Lagrange equations obtained from the Lagrangian 



L(x,y,a,x,y) 



-^Sx—a 



(-X' 



+ «')- r . 



K 



a— 3x 



p e 3hx + e 6x V 



V6 



(9) 



where we introduced the following dimensionless variable 

K 



y 



The equation dL/da = d(dL/da)/dt = leads to the constraint Eq.M). Fixing the 
gauge a = F(x,y), we can consider Eqs.(^|),(||) as the Euler-Lagrange equations obtained 
from the Lagrangian (0) under the zero-energy constraint (J|). 

In what follows we consider the potential of the form ([I]). As the system is 
symmetric under the transformation o — > —a, (p — > —ip, without a loss of generality we 
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will consider only the case a > 0. For the exponential potential ([I]) the Lagrangian (|9|) 
looks as follows 

L(x, y, a, x, y) = \e Zx ~ a (-x 2 + y 2 ) - ^e a - 3x [p e^ + V S X ~^] . (10) 
Dynamical system described by the Lagrangian of this form belongs to the class of 



pseudo-Euclidean Toda-like systems investigated in our previous papers |L2], [13], [14 



Methods for integrating of pseudo-Euclidean Toda-like systems are based on the 
Minkowski-like geometry for characteristic vectors («o + 3(/i — 1),0) G IR 2 and 
(a + 3, —6a) E IR 2 appearing when one puts the gauge a = a x with a =const. 
Here we do not describe the methods and refer to the above mentioned papers. 



3. General solutions 

A. a is arbitrary, h = 0. 

Here we suppose that the perfect fluid pressure p is equal to its density p (Zeldovich- 
type, or stiff matter). In this case the system with Lagrangian ( [TO]) is integrable for 
arbitrary parameter a in the so-called harmonic time gauge defined by 

a{x) = 3x. (11) 

Let us consider two different cases: o ^ 1 and a = 1. 
If a 1 we introduce the following variables 

u = —ax + y, v = x — cry. 

In the terms of u and v the equations of motion and the zero-energy constraint look as 
follows 

0, 



u 



V (l-a 2 )e 



6; 



• 2 i -2 

—U+V 



I -a 2 



Po + V e 



6ti 



We notice that for v we get the Liouville equation. Integrating this set of equations and 
inverting the linear transformation, we get the following general solution: 
the scale factor 

1/[3(<T 2 -1)] 



a = e 
the scalar field 



a 



f{t _ to)e eA(t-to) 



V6( 



111 



1/[3( CT 2 -1)] 



+ 2/o 



r( t_ to y(t-t ) 

where we introduced the function 

f{t - t ) = sinh(v/fi|t - t \)/VB, V (l - a 2 ) > 0, B > 0, 
= cosh(v / 5|t - t \)/y/B, Vo(l - (x 2 ) < 0, B > 0, 
B\t - t \)/VB, V (l - a 2 ) > 0, B < 0, 



(12) 
(13) 



sin 



If-fnl, V n (l 



a' 



> 0, B = 0. 
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The constant B is denned by 

B = A 2 + 3(1 - a 2 ) Po . 

The general solution contains 2 arbitrary constants t , A and 2 constants a ,y , obeying 
the constraint: 

e 6CT ^ = 3a 6 \V {l - a 2 )\. 

For the second case a = 1 the equations of motion and the zero-energy constraint 
with respect to the harmonic time gauge read in the old variables 

x = V Q S X ~ V \ 

y = Vo S x ~y\ 

-x 2 + y 2 = - l -[ P o + v,S*-y\ . 

We immediately find the integral of motion x — y — const = 2A. 

If Vo > the constant A is nonzero due to the zero-energy constraint. The general 
solution in this case for arbitrary Vq looks as follows 



a = e x = exp < (A + 



Po 
12A 



(t - t ) + V 



e 12A(t-t ) _ l 



+ X 



-A 



12A 



(t - t ) + V 



s 12A(t-t ) 



- 1 



Vo 



(14) 
(15) 



where the constants A, to are arbitrary and the constants Xo,yo obey the relation 
y - xq = 2At . 

For Vq < the additional solution appears (corresponding to A — 0, x — y) 



expjs^-to) - y(t-to) 2 + x } 



<P 



lna + -ln(-po/Vo) 
o 



where B and Xq are arbitrary constants. We remind that t is the harmonic time. It is 
connected with the cosmic time t c by the differential equation dt c = a 3 dt. 

B. a = 1 - h/2, 0<h<2. 
We notice that the model for a — 1/2 and h = 1 has been integrated in Here we 
study a more general case. We fix the time gauge choosing the following function a(x): 

a(x) = 3(1 — h)x. 

Then, with respect to the new variables u and v defined by 



u = exp 



3h 



(x-y) 



v = exp 



3h 



(x + y) 



the equations of motion and the zero-energy constraint may be written in a rather simple 
form 



u = 0, 



-h{2 - h)Vou^/\ 



uv 



Po 



+ Vou 2 ^' h ] 
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The set of equations is easily integrable. We obtain 

u = A(t — t ) > 0, 

3h 2 < 
V= 4A 2 

where A is an arbitrary nonzero constant and B is an arbitrary nonnegative constant, 
Vq has arbitrary sign. If Vq < 0, then, the following additional special solution arises 



p A{t 



to) + l^Vo[A(t - t )f" h ^ h + b\ > 0, 



u 



v = -h(2 — h)p 



h/(2-h)} 



T 2 -{t- t ) 



where the integration constant T / 0. Then, one easily gets the scale factor 



a = e 

and the scalar field 



[uv 



,l/(3fc) 



2/3 v 
In — . 



(16) 



(17) 



nh u 
C. ctV (1 - /i/2) 2 , 0</i<2. 
We introduce the following variables for the factorization of the potential in the 
Lagrangian (flI3|) 

tx = 3[<rx - (1 - V 2 )?/], 



V = 3[(V2-l)x + ay]-ln^|K,|/po- 

Then, the equations of motion and the zero-energy constraint for u and t> in the harmonic 
time gauge defined by Eq. flTTD look as follows 



u 



'-haA (e 2v + e) exp 



hau + (2 - h - 2a 2 )v' 



-Ac h(l - -)e 2v + (2-h- 2a 2 )e exp 



[l-h/2) 2 } (e 2 *+e) 



a 2 - (1 - h/2) 2 J ' 

( hau + (2 - h - 2a 2 )v ' 
\ a 2 - (l-h/2) 2 
\ hau + (2-h-2a 2 )v s 



exp 



a 2 - (1 - h/2f 



-ii 2 + v 2 + 3A [a 2 
where we denoted 

_ -2 — ti — -zo~ h(l—h/2) 

A = p 2[a ' 2 -ww\v \^* -d-v^i, 
e = sgn(Vfc). 

Let us express i> 2 from the zero-energy condition (|2"0"|) as follows 



(19) 
(20) 



-i 



3Ao[a 2 - (1 - /i/2) 2 ] (e 2 " + e) exp 



f/i(Jii + (2-/i-2cT 2 )x;' 
\ a 2 -(l-/i/2) 2 



Substituting ii, i; and w 2 into the relation 

d 2 u u — vdu/dv 
d 2 v v 2 



.(21) 



(22) 
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we obtain the following second order ordinary differential equation 



fu 
dv 2 





a 



h 2 /A) 



-,2c 



— e\ du 



<j- 



— + 



ha/2 



■ (23) 



-(l-h/2) 2 e 2v + ejdv a 2 - (1 - h/2) 2 

The procedure is valid if v ^ 0. 

The exceptional solution with v = appears only for the positive potential when 
a 2 > 1 — h/2 and < h < 2. In the terms of the cosmic time t c it reads 

* 2/[3(2-h)] 

-a(2 - h)\tr - t°J I 



A(a 2 - (1 - h/2) 




(2-h)V a\t c -t° c \ 



(24) 
(25) 



It should be mentioned, that the set of the equations ([18| ) - (|20|) does not admit 
static solutions ii = v = as well as the solutions with u = dbu. So, using the relations 
(f2~l"l) and (|22|) we do not lose any solutions except, possibly, the exceptional solution 



Let us suppose that one is able to obtain the general solution to the equation (^3|) in 
the parametrical form v = v (r), u = u(r), where r is a parameter. Then, we obtain the 
scale factor a = exp[x] and the scalar field cp = (V&/ n)y as functions of the parameter r. 
The relation between the parameter r and the harmonic time t may be always derived by 
integration of the zero-energy constraint written in the form of the following separable 
equation 



dt 2 



dr J 



dr J 



3A [a 2 - (1 - h/2) 2 } (e 2v + e 



■ exp 



hau + (2-h-2a 2 
a 2 - (l-h/2) 2 



dr 2 



(26) 



Thus, the problem of the integrability by quadrature of the model is reduced to 
the integrability of the equation (p3|). For du/dv it represents the first-order nonlinear 
ordinary differential equation. Its right hand side is the third-order polynomial (with 
coefficients depending on v) with respect to the du/dv. An equation of such type is 
called Abel's equation (see, for instance, |T5|],||T6|). 

First of all let us notice that the equation (|23|) has the partial solution u = ±w+const 



that make the relation (^T|) singular. However, as was already mentioned, the set 
of Eqs. (|l8D -(E2~0D does not admit the solutions with it = ±v. Existence of this 
partial solution to the Abel equation (^) allows one to find the following nontrivial 
transformation 



'2i> 



XdY 



— e 



Y dX 



u 



v + In 



Y 
X 



+ lnC 



6 = ±1, C>0, 



that reduces it to the so-called generalized Emden-Fowler equation 



d 2 Y 
dX 2 



sgnfcr 2 - (1 - V2) 2 ] s 



dY_ 

dX 



Yin j^r 



(27) 
(28) 

(29) 



Friedmann universe with perfect fluid and scalar field 



where the constant parameters /, m and n read 



I 



2(5a - 1 + h/A) 



5a + 1 - h 



m 



n 



-m — 3. 



(30) 



5a-l + h/2 ' "" 5a + l-h/2 : 

In the special case I = Eq. (p9|) is known as the Emden- Fowler equation. 

There are no methods for integrating of the generalized Emden-Fowler equation 
with arbitrary independent parameters I, m and n. However, the discrete-group methods 
developed in [Tj| allow to integrate by quadrature 3 two-parametrical classes, 11 one- 
parametrical classes and about 90 separated points in the parametrical space (I, m, n) 
of the generalized Emden-Fowler equation. Further we consider the following integrable 
classes. 

CI. a = l- h/A, 0<h<2, 5=1. 
The parameters I and m given by Eq. (|30"D are the following 

2h 



I = 0, m 



-1 + 



■1,1]. 



(31) 



8-3/i 

The general solution to the generalized Emden-Fowler equation (|29"D with these 
parameters reads 

t _ „ 1 



Y 



Fjr) 



>0, X 



F m (r) 

where we introduced the following function 

2 



> 0, 



Fm.(r) 



± 



-T 



m+l 



-1/2 



dr + C 2 . 



(32) 



.m + l 

The variable r changes on the interval which follows from G m (r) > 0, where we used 
the function 

" F m (r) 



G m (r) 



tFL(t) 



equal to the right hand side of Eqs.(|2"T|) with substitutions of X,Y and 
using Eqs. ([27l) ,(|2~8D we find the scale factor 

4(4-h) J_ 

a = a r 3h (s-3h) G%{t) 



and the scalar field 



<P 



hi 



8(2-h) _2_ 
T 3M8-3h) G%(t) 



+ 2/o 



(33) 
Finally, 

(34) 
(35) 



where 



4(4-h) 4,(2- h) 

a Q = C , 3"(8-3h)(|^ |/p )3M8-3h) } y Q 



In 



8(2-h) 2(4-h) 
(J3h(8-3h) (| Vo|/po) 3/l(8 " 3,l) 



The relation between the variable r and the cosmic time t c is the following 



4(2-fe)(4-3h) 2(2-h) 

dt 2 c = JorWft, - (r) \F' m {r)\ 2 dr 2 , 



(36) 



where 



16p 



T^y,..,^) 2 „4(2- h )(4-h) 



y |Kwk)CT^r/[3/j(8 - 3/i]. 
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The general solution (|34|) , (|35|) (with I and m given byfl3~l|)) contains 3 arbitrary constants 
C, C\ and C 2 as required. 

C2. a = |1 -h\, < h < 2, h=£ 1,4/3,5 = sgn(l - h). 
Now we have the following parameters in the generalized Emden-Fowler equation 

h 



I 



m 



-1 + 



4-3/i 

Its general solution looks as follows 

t ..1 



X 



>0, Y 



Rm\r) Rm{r) 
where the function R m (j) is defined by 

-2e 



>0, 



Rm{r) = ± 



± 



-m-2 



\m + 2\ 
elnT 2 + d 



-1/2 



-1/2 



dr + C 2 , m 



dr + C 2 , m ^ -2, 
-2. 



The variable r changes on an interval, where 



S m (r) 



Rm(r) 
tRL(t) 



> 0. 



Finally, we get 



4(l-h) J_ 

a = a T 3h ^-3h)S^(j) 



where 



<P 



a 



-^sgn(l — h) (in 

K I 



2(2-h) 1 



4(l-h) 2(2-fc) 
C3h(3h-4) (|V r |/p )3h(3h-4) ( = 111 



+ 1/0 



2(2-h) 4(l-h) 
(73fe(3fe-4) (lyj /p Q )3h(3h-4) 



The variable r and the cosmic time t c are related by 



dt 



4(2-3ft) 



/7 t 3,1 S 



where 



(2-fc) z 



4(2-h)(l-h) 



(37) 



(38) 
(39) 

(40) 

(41) 
(42) 



(43) 



U = 4p h{:ih - i) \V \ t (^C '"(3/1-4) 7(3/i|3/i-4|). 

So, this general solution is given in the parametrical form by (fill) , (f42|) with / and m 
from (|37|). The transition to the cosmic time may be done by solving ([43|). 



4. Properties of solutions 

Here we study properties of the obtained exact solutions for positive potentials, though 
our solutions are valid for any sign of the potential. As the system is symmetrical under 
the time reflection t — > — t, without loss of generality we only consider expanding near 
the singularity cosmologies with the Hubble parameter H = xexp(— a) > 0. 
We introduce the following notation for the relative energy densities 

*V K 2 (d^/dt c ) 2 _ #V(<p) 

p ~ 3# 2 ' vK ~ 6# 2 ' VP ~ 3# 2 • 
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Due to the constraint (|j) the relative energy densities obey the relation 
Also we introduce the scalar field barotropic parameter 

Al. < a < 1, h = 0. 
The general solution is given by Eqs.(p|),(|13|) for Vo(l — o -2 ) > 0, B > 0. Near the 
initial singularity (t c — > +0) the scale factor is in the main order a oc t l J 3 . There exists 
the special solution for A = —o\f~B with Vl p — > 1 as t c — > +0. It corresponds to the 
fluid-dominated solution mentioned in 0, |9|, [10]. It gives p v / p,p — > 1 (stiff matter) as 
£ c -> +0. 

For all remaining solutions we obtain Q p — > 1 — — > tylx, Q v p — > as 
t c — > +0, where = [o\[B — A) 2 /(VB — a A) 2 . Then the barotropic parameter p^/ p v 
tends to —1, i.e the scalar field is vacuum-like (de Sitter) near the initial singularity. 

There is an unique late-time attractor in the form of the scalar field dominated 
solution. It corresponds to B — 0, and p = in formulas fll2|),(|L3|) presenting the 
general solution. The attractor may be written down for the cosmic time t c 



a = a Q 



t l/(3<x 2 ) 



2/3 



In t c + In 



'3V 
1^ 



It is easy to see that for this solution 



a, = o, n 



a 



<pp 



l — a 



Pv 



1o l - 1. 



This attracting solution according to |L0[ may be called kinetic-potential scaling. For 
o 2 < 1/3 all solutions provide the power law inflation at late times. 
A2. l<a,h = 0. 

The general solution is given by fll~2|) , (|13|) for Vo(l — u 2 ) < 0, B > 0. Behavior near 
the initial singularity is the same for both general and special solutions with the only 
difference that the constant Vt Q K is the following Vt° K = (ay/B + A) 2 /{\fB + a A) 2 . 
At the final stage of evolution as t c — > +oo we obtain Q p — > 1 — ^lL K , Q v k 



a oc t^ 3 , where 



(aVB - A) 2 /(VB -a A) 2 . Such 



behaviour of the scalar field, when it adopts the usual perfect fluid equation of state and 
its energy- density scales with that of the perfect fluid, is called scaling (or sometimes 
" tracking" ) . 

B and C2 in regions la and lb. 
Behaviour near the initial singularity was described in |7|, |], |10j using qualitative 
methods. There exists the fluid dominated solution with a oc ^/[3(2-/i)] anc j o ]_ 



Friedmann universe with perfect fluid and scalar field 



12 



£LS t r 



-0. All remaining solutions describe the domination of the kinetic contribution 



of the scalar field: fL 



u p - 0, Q^k -> 1, Q vP -> 0, p v /p v — 
The behaviour at late times is the same as one for Al. 



a oc 



tV\ 



CI and C2 in regions Ila and lib. 
Behaviour near the initial singularity is the same as in the previous case. The late-time 
attractor for h e (0, h) is the special solution described by Eqs. (|2"3|) and (|2"5|). For these 
solution we have 

2-h „ {2-hf h(2-h) p v 

liiaP — : — ~ , — — 



4(T 2 > "VP 4a 2 

This is a typical scaling behaviour. 

One of the examples of this behaviour may be given explicitly: the exact solution 
of the class CI with h = 2 and a = 1/2 for the positive potential reads with respect to 
the cosmic time 

1/3 
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1 + A, 
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where A = po is the cosmo logical constant. The solution contains 3 integration 
constants: arbitrary £° positive a and A obeying \A\ < 1. The late-time attractor 



corresponds to A = 0. For this attracting solution we get: Q p — > 1, 
0, Pip/ptp— 1, a (x t x J 3 exp J A/3t c , H — > J A/3 as t c — >■ +oo. 
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